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Experimental  investigations  and  numerical  modeling  of  optical 
phase  modulation  by  dual-frequency  nematic  liquid  crystals 

Mikhail  I.  Barnik,  Vladimir  V.  Lazarev,  Serguei  P.  Palto 

Shubnikov  Institute  of  Crystallography  of  Russian  Academy  of  Sciences,  Leninsky  pr.59, 

119333  Moscow,  Russia 

Experimental  study  and  numerical  simulations  of  optical  phase  modulators  based  on 
nematic  liquid  crystals  characterized  by  frequency  sign  inversion  of  the  dielectric 
anisotropy  has  been  performed.  The  received  results  point  out  an  extreme  role  of  back¬ 
flow  effects  and  electro-hydrodynamic  instability,  which  make  restrictions  for  creating 
high-speed  phase  modulators. 

1.  Introduction 

For  a  long  time  nematic  liquid  crystals  (NLCs)  are  used  in  adaptive  optics  for  wavefront 
shaping  and  control  (see  [1-14],  for  instance).  The  optical  phase  modulation  realized  in  such 
devices  is  due  to  optical  retardation,  which  appears  as  a  result  of  field-driven  NLC  director 
(optical  axis)  reorientation.  For  instance,  the  planar  LC  director  field  distribution  characterized 
by  in-plane  orientation  of  the  optical  axis  can  be  switched  by  electric  field  to  the  homeotropic 
director  field  distribution,  when  the  optical  axis  becomes  normal  to  the  LC  layer.  In  case  of 
positive  dielectric  anisotropy  of  the  LC  (sa=s\\—s±>0,  where  £\\  and  s±  are  the  dielectric 
permittivity  along  and  perpendicular  to  the  LC  director  respectively)  the  latter  is  possible  by 
using  an  electric  field  applied  along  the  LC  layer  normal.  In  this  example  the  switching  speed  to 
the  homeotropic  state  can  be  raised  by  increasing  the  electric  field  strength.  For  typical  LC 
parameters  and  reasonable  strength  of  the  electric  field  of  about  10  V/pm  the  field-on  switching 
time  can  be  as  low  as  lOOps.  However,  the  back  relaxation  to  the  ground  planar  state  is  still 
defined  by  viscous-elastic  properties  of  LC  material,  so  the  relaxation  speed  is  restricted.  Even 
for  low  viscosity  materials  (with  the  rotational  viscosity  y~0.05  Pa  s,  and  elastic  coefficients 
Ki,2,3~10pN)  at  LC  layer  thickness  of  5pm  the  relaxation  time  is  of  about  15ms.  Quite  similar 
situation,  but  with  opposite  direction  of  the  optical  axis  driving,  takes  place  if  the  initial  LC 
director  field  distribution  is  homeotropic,  but  the  LC  has  the  negative  dielectric  anisotropy  (sa< 0). 
Thus,  one  of  the  main  drawbacks  of  the  LC-based  modulators  that  make  restrictions  for  the 
modulation  speed,  is  due  to  existence  of  the  free  viscous-elastic  relaxation.  Fortunately,  there  are 
liquid  crystal  materials,  which  allow  excluding  the  free  relaxation  stages  from  the  modulation 
process.  These  are  so-called  dual-frequency  liquid  crystals  (2f-LCM)  [15-20],  which  show 


inversion  of  the  dielectric  anisotropy  sign  at  some  frequency,  which  is  typically  located  in  sound- 
frequency  band.  Using  the  dual-frequency  driving  one  can  realize  field-driven  reorientation  of 
the  LC  optical  axis  to  both  homeotropic  and  planar  states.  There  are  many  examples  of  2f-LCMs 
applications  for  the  optical  phase  modulation  [21-26],  Due  to  possibility  of  excluding  the  free 
relaxation  stages  from  the  modulation  process  it  is  believed  that  using  the  2f-LCMs  will  allow 
creating  very  fast  phase  modulators  [27,28], 

Recently  [29]  we  have  studied  the  dynamical  properties  of  2f-LCMs  phase  modulators. 
The  stable  switching  to  any  desired  state  characterized  by  phase  retardation  within  the  range 
from  0  to  2k  has  been  achieved  for  times  less  than  1  ms.  This  switching  was  realized  for  layers 
of  thicknesses  larger  than  7  pm.  For  example,  to  obtain  such  phase  modulation  parameters  for  a 
layer  of  the  2f-LCM  as  thick  as  8  pm,  it  is  sufficient  to  use  the  low-frequency(2  kHz)  driving 
voltage  of  50  V  and  high-frequency  voltage  (35  kHz)  of  30  V.  We  have  also  assumed  that  high- 
frequency  driving  induces  extra  back-flow  coupling,  which  then  provokes  experimentally 
observed  hydrodynamic  instabilities  restricting  high-speed  performance  of  the  modulators. 

In  this  paper  we  clarify  reasons  for  restrictions  of  high-speed  modulation  in  more  details. 
In  addition  to  the  experimental  study  we  show  the  results  of  the  numerical  simulations  of  the  2f- 
LCMs  switching.  These  numerical  simulations  agree  with  the  experimental  observations  and 
give  explicit  confirmation  for  the  principal  role  of  back-flow  hydrodynamic  coupling,  which 
makes  restriction  in  creating  high-speed  optical  phase  modulators.  The  simulation  results  point 
out  the  principal  LC  parameters  to  be  improved  in  order  to  get  highest  performance  of  the  2f- 
LCM  phase  modulators. 

2.  Experiment  and  Results 

Experimental  data  were  received  on  a  specially  developed  2f-LCM  mixture  called  AP- 
99-1.  This  mixture  is  characterized  by  following  physical  parameters  at  room  temperature  (T  = 
21°C):  sa  (1  kHz)  =  2.4,  sa  (40  kHz)  =  -2.6,  f  =  9  kHz,  optical  anisotropy  An  =  «//  -  n±  =  0.25  («// 
=  ne  and  n±  =nQ  are  refractive  indices  along  and  perpendicular  to  the  LC  director,  respectively), 
splay  and  bend  elastic  constant  are  equal  K\  =  9.7  pN  and  K]  =  16.5  pN,  respectively.  The 
temperatures  of  phase  transitions  are  found  to  be  as  follows:  from  nematic  to  smectic  TNs=<-7°C 
and  from  isotropic  to  nematic  Ti_n=74°C  at  cooling  and  from  nematic  to  isotropic  TN-i  =78°C  at 
heating. 

The  electro-optical  investigations  were  made  on  cells  with  planarly  aligned  LC  layers. 
The  LC  alignment  was  due  to  the  boundary  conditions  providing  the  director  pretilt  angles  of  the 
same  sign  at  opposite  surfaces.  These  boundary  conditions  were  controlled  by  rubbing  of  the 


polyimide  alignment  layers  in  opposite  ("antiparallel")  directions  that  gives  a  pretilt  angle  of  |0| 
~3°-4°  [30].  The  cell  of  such  a  director  arrangement  is  referred  below  as  a  cell  of  homogeneous 
alignment,  or  "0-aligned"  cell.  For  the  0-aligned  cells  the  director  field  is  uniform  accrosst  the 
thickness  of  LC  layer.  The  cells  were  placed  between  crossed  polarizers.  The  He-Ne  laser  was 
used  as  a  light  source.  The  dual-frequency  modulated  pulse  packets  consisting  of  a  sequence  of 
low  frequency  (1  kHz)  square- waveform  pulses  and  high  frequency  (35  kHz)  sinusoidal  waves 
were  used  for  driving.  The  details  of  the  method  for  measuring  the  optical  phase  modulation 
properties  were  described  in  [29], 

Two  geometries  distinguished  by  LC  director  orientation  (i.e.  rubbing  direction  R  of 
alignment  layers)  with  respect  to  the  direction  P  defining  the  axis  of  the  input  polarizer  were 
used  in  our  experiments:  i)  “RP45”,  when  the  rubbing  direction  turned  at  an  angle  of  45°; 
(ii)“RPo”,  when  the  rubbing  direction  is  at  an  angle  of  0°.  In  all  the  measurements  the  input  and 
output  polarizers  were  crossed.  For  RP45  geometry,  in  case  of  non- twisted  LC  director  field 
distribution,  the  intensity  /  of  the  outgoing  light  and  the  phase  retardation  AO  are  related  as: 


I  =  I0  sin2 


AO^| 

{Tj 


(1) 


where  Io  is  approximately  (neglecting  reflections  and  non-ideality  of  the  polarizers)  half  of  the 
intensity  of  the  incoming  nonpolarized  light  and 


AO  =  ^-dAn, 
X 


(2) 


where  X  is  wavelength  and  d  is  LC  layer  thickness.  Note  right  here  that  for  RPo  geometry  the 
electro-optical  response  can  be  non-zero  only  if  a  twist  deformation  appears  in  the  director  field 
distribution. 

Fig.  1  demonstrates  the  electro-optical  response  for  two  different  cells  with  LC  layer 
thickness  of  7.6  pm  and  10  pm.  For  these  thicknesses  values  the  field-off  states  correspond  to 
phase  retardations  of  about  671  and  871,  respectively.  Maximum  values  of  transmission  (T=l)  are 
achieved  at  odd  multiple  of  71  phase  retardations.  The  driving  waveform  is  shown  at  the  bottom 
in  the  Fig.  1 .  Both  low  and  high  frequency  driving  waveforms  are  of  1  ms  in  duration. 


a) 


b) 


Fig.  1.  Electro-optical  response  measured  for  the  0-aligned  AP-99-1  layers  in  RP45  and  RPo 
geometries  at  Ulf  =  50V  and  T  =  22  C):  a)  d  =  7.6  pm,  curve  1-4  are  for  Uhf  =22,  30,  40  and  0 
V,  respectively;  b)  d  =10  pm,  curve  1-3  are  for  Uhf  =  30,  38  and  51V,  respectively.  The  lower 
trace  represents  the  applied  voltage  waveform. 


As  it  is  seen  in  Fig.  1,  the  low  frequency  driving  waveform  of  about  50  V  in  amplitude  is 
sufficient  to  achieve  almost  homeotropic  director  field  distribution  that  provides  the  phase 
retardation  AO  to  be  very  close  to  zero.  In  the  dynamic  mode  the  phase  retardation  of  2tc  not 
always  can  be  achieved,  despite  the  two  layers  are  capable  to  provide  a  phase  retardation  of  671 
and  871  respectively  in  the  quasi-static  regime.  The  latter  is  proved  by  absence  of  100%  intensity 
modulation  in  the  electrooptical  response.  The  amplitudes  of  high  frequency  voltage  Uhf  are  30 
and  38  V  for  LC  layer  thickness  of  7.6  pm  and  10  pm  respectively.  At  first  glance  with 
increasing  the  high  frequency  amplitude  Uhf  at  a  fixed  low  frequency  voltage  Ulf  the  LC 
director  field  distribution  after  fixed  time  interval  should  become  closer  to  the  planar  state,  and 
the  phase  retardation  have  to  be  increased.  However,  one  can  see  that  after  exceeding  some 
certain  value  of  Uhf,  the  maximum  transmission  in  RP45  geometry  begins  to  fall  down  and  100% 
intensity  modulation  is  not  reached.  At  the  same  instant,  the  step-wise  growth  of  the  transmission 
appears  in  RP0  geometry.  A  non-zero  transmission  in  RP0  geometry  can  only  be  explained  if  LC 
director  field  is  subjected  to  the  twist  deformation.  By  additional  experiments  we  have 
established  that  the  observed  behavior  takes  place  just  before  the  electrohydrodynamic  (EHD) 
instability,  which  appears  at  larger  high-frequency  fields.  The  EHD  is  a  result  of  complex 


hydrodynamics  developed  in  LC  layer,  when  high  frequency  driving  voltage  exceeds  a  certain 
threshold  value.  At  voltages  higher  than  threshold  values  the  EHD  instability  is  characterized  by 
appearance  of  periodic  distortions  of  LC  texture.  In  RPo  geometries  the  EHD  instability  looks 
like  striped  (black-bright)  domains  oriented  along  the  rubbing  direction.  At  further  growth  of 
Uhf,  regular  hydrodynamic  flows  become  turbulent  and  accompanied  by  strong  light  scattering 
[29], 

The  experimental  study  was  performed  for  LC  layers  of  thickness  varied  in  a  range  of 
2.5-12.5  pm.  In  Fig.  2  we  summarize  the  experimental  results  on  thickness  dependencies  of  the 
high-frequency  driving  voltage  Uhf  for:  i)  EHD  instability  appearance  (curve  1);  for  stable  271- 
modulation  (curve  2);  iii)  for  stable  7i-modulation  (curve  3).  The  low  frequency  voltage  Ulf  in 
all  these  experiments  is  maintained  at  a  constant  value  of  50  V.  The  plot  area  above  the  curve  1 
corresponds  to  EHD  instability.  Due  to  EHD  instability  restrictions  the  fast  (500  Hz)  stable  2tz 
phase  modulation  can  only  be  realized  if  LC  layer  thickness  is  larger  than  7  pm.  For  thinner 
layers  within  a  thickness  range  from  3  pm  to  6  pm  a  maximum  phase  modulation  of  n  can  only 
be  achieved  (see  Fig.  2,  curve  3). 


Fig.  2.  Thickness  dependence  of  high  frequency  (35  kHz)  voltage  Uhf  at  which  a  dynamic  phase 
retardation  change  of  n  (3)  and  2n  (2)  is  still  achieved  (100%  light  modulation  efficiency  is 
observed).  The  curve  1  shows  dependence  of  the  threshold  voltage  of  dynamic  phase  modulation 
breakdown.  The  data  are  for  the  0-aligned  AP-99-1  layers  at  T  =  22  C 


Fig.  3.  A  sketch  of  the  double-layer  liquid  crystal  cell  and  its  connection  to  the  net  source.  The 

arrows  show  the  directions  of  rubbing. 

Double-layer  LC  cells  can  be  used  to  expand  the  dynamic  range  of  fast  optical  phase 
modulation.  In  this  case  the  same  speed  at  equivalent  modulation  is  achieved  at  lower  values  of 
the  high-frequency  voltage,  when  LC  is  not  subjected  to  EHD  instability.  A  double-layer  cell 
design  and  its  electric  connection  to  the  net  source  are  schematically  shown  in  Fig.3.  The  cell 
assembly  consists  of  three  stacked  glass  plates  with  two  gaps  between  them,  which  are  filled  by 
AP-99-1  LC  material.  The  inner  surfaces  of  the  two  outer  plates  and  the  both  surfaces  of  the 
middle  plate  have  ITO  electrodes,  coated  with  alignment  polyimide  layers,  which  are  rubbed  in 
the  directions  designated  by  arrows  in  the  Fig.  3  to  provide  planar  homogeneous  LC  alignment. 

Fig.  4  shows  oscillogramms  of  electro-optic  response  for  the  double-layer  cells  with 
different  LC  layers  thickness  of  6.8  and  7.3  pm,  when  the  electric  circuit  providing  the  parallel 
connection  of  layers  to  the  voltage  source  is  used  (see  Fig.  3).  Thus,  each  layer  is  driven  by  the 
same  periodic  pulse  sequence  consisting  of  the  low  frequency  (1  kHz)  square-waveform  at  fixed 
voltage  values  Ulf  =  50  V,  and  the  high  frequency  35  kHz  sinusoidal  waveform  at  different 
values  Uhf-  Firstly,  it  is  necessary  to  note  that  maximum  transmission  of  the  double-layer  cell  in 
case  of  RP45  geometry  is  approximately  twice  less  than  that  of  the  single-layer  cell  (compare 
with  Fig.  1).  This  is  an  effect  of  increasing  from  2  to  4  in  the  number  of  ITO  electrodes,  which 
absorb  and  reflect  some  portions  of  light  (transmission  coefficient  of  a  single  plate  with  ITO 
electrode  we  used  is  as  large  as  0.8).  The  problem  of  reduced  transmission  can  be  resolved  to  a 
great  extent  by  known  technological  ways:  using  a  less  absorbing  electrode  coating, 
antireflective  coverings  and  so  forth.  Another  feature  of  the  double-layer  cell  operation  is  an 
increased  light  leakage  in  the  dark  state,  when  the  LC  director  is  driven  to  the  homeotropic  state. 
The  latter  is  due  to  increased  number  of  LC  interfacial  layers,  which  are  never  reoriented  to  the 
ideally  homeotropic  states  because  of  the  strong  surface  anchoring.  This  light  leakage  in  the  dark 
state  can  be  considerably  reduced  by  increasing  the  voltage  amplitude  of  low  frequency  driving 


and/or  by  weakening  the  anchoring  conditions  at  the  interfaces.  Another  way  is  using  of  external 
retarder  fdms  to  provide  optical  compensation  of  these  quasi-homeotropic  interfacial  LC  layers. 
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Fig.  4.  Electro-optical  response  of  the  0-aligned  AP-99-1  layers  in  a  double  layer  (6.8  and  7.3 
pm)  cell  measured  in  RP45  and  RPo  geometries  at  Ulf  =  50V  (T  =  22  C):  curve  1-4  are  for  Uhf 
=0,  22,  28  and  32  V,  respectively 


The  basic  positive  result  obtained  on  the  double-layer  cell  is  the  following.  It  is  clearly 
seen  from  the  Fig.  4  that  the  equal  in  value,  say  2k,  phase  retardation  change  for  the  double-layer 
cells  and  ordinary  single-layer  cells  is  reached  at  different  high  frequency  voltages.  The  voltage 
is  by  a  factor  of  -  1.4  less  for  the  double-layer  design  compared  to  single-layer  one.  In  particular 
case  the  high  frequency  voltage  magnitudes  are  22  V  for  the  double-layer  design  and  -32  V  for 
the  single-layer  cell.  Owing  to  such  a  marked  reduction  of  driving  voltage,  it  is  possible  to 
realize  stable  switching  of  phase  retardation  even  in  the  range  from  0  to  3n,  while  for  the  single¬ 
layer  cell,  at  the  same  voltages,  the  undistorted  fast  phase  retardation  switching  is  hardly 
realized  in  the  dynamic  range  from  0  to  2n. 


3.  Numerical  simulation 


We  made  comprehensive  numerical  analysis  of  the  light  polarization  state  dynamics  for 
the  two  virtual  systems  to  clarify  reasons  for  restrictions  of  speed  and  value  of  phase  modulators 
in  case  of  using  the  dual-frequency  LC  materials.  For  the  first  type  of  the  system  the 
hydrodynamic  coupling  was  suppressed,  i.e.  the  dynamics  of  the  LC  director  field  deformation 
has  been  calculated  only  on  account  of  the  rotational  viscosity,  while  in  the  second  case  the  set 
of  dynamic  equations  is  extended  to  take  into  account  the  hydrodynamics.  In  the  latter  case  the 
viscous  properties  are  described  by  five  independent  Leslie  coefficients. 

The  simulations  have  been  performed  using  LCD  TDK  2.0  software  package,  which  was 
created  by  one  of  participants  (P.S.P.)  of  this  project  for  solving  different  problems  in  field  of 
electro-optical  properties  of  LC  devices. 

3.1.  Simulated  Set-Up 

The  scheme  of  the  simulated  set-up  is  shown  in  Fig  .5.  The  LC  cell  is  placed  between  two 
crossed  polarizers  with  the  transmission  axes  oriented  at  an  angle  of  45  degrees  with  respect  to 
the  x  axis  of  the  laboratory  frame  xyz.  The  polarizers  are  assumed  to  be  of  high  optical 
performance.  In  calculations  the  polarizers  are  simulated  by  plates  of  thickness  of  2  pm  made  of 
material  providing  the  dichroic  ratio  of  k|/k j  =100  (k||  =3  pm'1,  k j  =0.03  pm'1). 


Fig.  5.  Scheme  of  simulated  electro-optical  set-up. 


The  LC  layers  of  thickness  of  5  jam  are  simulated,  and  they  are  assumed  to  be 
homogeneously  aligned,  so  at  the  layer  boundaries  the  LC  director  has  almost  planar  orientation 
along  the  x  axis  (small  pretilt  angle  of  4  degrees  typically  realized  in  practice  is  defined  in 
calculations).  To  study  the  stability  of  the  director  field  with  respect  to  the  hydrodynamics  we 
have  also  introduced  small  azimuthal  angular  difference  of  one  degree  in  easy  axes  orientations 
at  the  LC  layer  boundaries.  This  difference  typically  exists  in  experimental  cells  due  to 
technological  inaccuracies  (rubbing  and  cell  assembling  procedure).  Some  of  LC  parameters  are 
fixed  in  simulations,  while  the  others  are  variable.  Around  the  fixed  quantities  are  the  elastic 
constants  (Ki=\ 0  pN,  K2= 5  pN,  K3=\ 5  pN),  rotational  viscosity  y=  0.5  Pa  s,  optical  anisotropy 
An  =  0.2,  low-frequency  dielectric  anisotropy  ea  =  +5,  high-frequency  dielectric  anisotropy  sa  =  - 
4.  The  used  values  are  close  to  ones  of  typical  dual-  frequency  materials  used  in  experimental 
part  of  this  work,  but  they  do  not  match  them  exactly.  The  latter  is  not  of  principal  significance 
for  the  general  features  of  the  discussed  problem. 

3.2  Basic  equations  for  calculations  of  the  dynamics  of  the  director  field  under 
external  electric  field 

We  are  considering  an  LC  layer,  which  is  homogeneously  aligned  in  the  xy  plane 
directions  of  the  laboratory  frame  system.  The  external  electric  field  can  be  applied  in  the  z- 
direction  that  is  perpendicular  to  the  LC  layer  plane.  Thus  the  inhomogeneous  director  field 
distribution  takes  place  only  along  the  z-direction,  and  the  mathematical  problem  is  one¬ 
dimensional. 

3.2.1.  Equations  in  case  of  neglecting  the  hydrodynamic  coupling 

In  cases  we  neglect  the  hydrodynamic  coupling  the  Euler-Lagrange  formalism  results  in 
the  following  set  of  equations  describing  the  dynamics  of  the  LC  director  field  distribution 
n(z,  t )  =  {nx  (z,  t),  ny  (z,  t),  nz  (z,  t)) : 


dn  dF  d  dF 

—  =  Anx - + - , 

dt  dnx  dz  dnx  ' 

dn  ,  dF  d  dF 

— -  =  A  nv - + - 

dt  y  dny  dz  dny  1 

dn  dF  d  dF 

— -  =  A n7 - + - , 

dt  dnz  dz  dnz ' 


(3) 

(4) 


r 


(5) 


'  d 

where  yis  the  rotational  viscosity,  X  is  Lagrange  multiplier  and  nx  =  — n  .  Typical  values 

dz 


of  the  rotational  viscosity  for  the  dual-frequency  materials  are  in  a  range  of  0.4  -  0.7  Pa  s.  In  the 
case  of  nonchiral  LCs  the  free  energy  density  F  is  expressed  in  terms  of  both  elastic  and  electric 
field  energy  densities  as  follows: 


F  =  '2+K2 i(nynx  ~nxny  ')2  +Ki3[nz  2(ny  '2 +nx  ,2)+(nxnx  '+nyny  ')2| 
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where  Kti  are  the  elastic  modules  responsible  for  different  types  of  the  deformations.  The  electric 
field  contribution  is  expressed  in  terms  of  electric  displacement  Dz,  which  remains  independent 
of  z  due  to  negligible  conductivity  of  the  LC  layer.  The  relationship  between  the  electric 
displacement  and  applied  voltage  V  is: 


V  =  |  E(z)dz  =  | 


dz , 


£l  ®(1  +£^£^nz2) 


(7) 


For  dual-frequency  materials  the  dielectric  anisotropy  sa  is  frequency  dependent.  It  is  positive  at 
low-frequencies  (typically  less  than  a  few  kHz)  and  negative  if  the  high-frequency  field  is 
applied  ( typically  f>  5-10  kHz) . 

3.2.2.  Equations  in  case  of  the  hydrodynamic  coupling  is  taken  into  account 


On  account  of  Leslie  theory  [31]  the  total  set  of  equations  becomes 
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dnx  .  dF  d  dF 

y  — -  =  Anr - + - ®2nzK  > 

dt  onv  dz  onv 


dn  ,  dF  d  dF 

y  — -  =  A  nv - + - 

dt  y  dny  dz  dny ' 


a2nzVy 


dn  dF  d  dF  .  . 
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dt  dnz  dz  dnz ' 

where  VXti.  are  the  velocity  components  of  the  hydrodynamic  flow  and  a,  -  the  viscosity 
coefficients,  which  are  also  known  as  Leslie  coefficients.  The  rotational  viscosity  is  also 
expressed  in  terms  of  Leslie  coefficients  as  y=a^-a2.  According  to  theoretical  and  experimental 
studies  the  values  of  a\  and  <23  are  much  less  than  the  values  of  the  other  coefficients,  so  they 
are  neglected  in  current  simulations. 


3.3.  Basics  of  the  approach  for  the  optical  calculations 

The  optical  calculations  are  based  on  exact  solution  of  Maxwell  equations  for  one- 
dimensionally  inhomogeneous  anisotropic  media.  The  Berreman  4x4  matrix  approach  [32]  with 
numerical  algorithm  described  in  [33]  is  used  to  calculate  the  transmission,  reflection  and 
corresponding  light  polarization  states  for  the  system  shown  in  Fig.  5.  In  case  of  the  light 
polarization  state  calculations  the  second  (output)  polarizer  is  not  taken  into  account. 

If  the  media  is  homogeneous  in  the  xy-plane  the  system  of  six  Maxwell  equations  is 
reduced  to  the  system  of  only  four  independent  equations,  and  the  wave  propagating  equations 
can  be  written  in  terms  of  the  first  order  differential  equations  as 


-T  =  -A'T,  (13) 

dzc 


where  the  components  of  matrix  A  are  expressed  in  terms  of  dielectric  tensor  components  [32, 
33],  and  the  electromagnetic  field  column  T  is  defined  through  the  electric  (Exy)  and  magnetic 
field  (Hxy)  components  as  follows: 

f  £  ^ 


¥  = 


H„ 


k-HxJ 


(14) 


It  is  evident  that  the  solution  of  (13)  for  an  homogeneous  optically  anisotropic  slab  of 


thickness  h  is 


¥(/*)  =  exp (ia>  h\  /  c)'F(O)  =  P(h)W(0) , 


(15) 


where  'F(O)  and  'F(Zj)  are  the  electromagnetic  field  columns  at  the  input  (z= 0)  and  output  (z=h  ) 
of  the  slab  respectively.  Actually,  the  problem  is  reduced  to  calculation  of  the  exponent  of  the 
matrix  A  and  finding  the  matrix  P(h).  In  case  of  multilayered  optical  systems  one  needs  to 
calculate  the  product  of  matrices  Pi(hi)  of  the  individual  layers  the  system  consists  of.  Similarly, 
in  the  case  of  the  inhomogeneous  LC  layer  the  problem  is  solved  by  space  discretization 
(sampling)  of  the  whole  layer  thickness  along  the  z-axis  and  finding  the  product  of  the 
corresponding  matrices.  The  details  on  how  to  efficiently  evaluate  the  field,  transmission  and 
reflection  are  described  in  [33], 

The  analysis  of  the  polarization  state  is  performed  in  a  frame  of  the  Stokes  vector  concept 
using  Poincare  sphere  presentation  after  the  output  field  is  calculated.  By  expressing  the  output 
electric  field  components  in  form 


Ex  =  Ax  cos(&tf  -kz  +  8 x ) 

Ey  =  Ay  cos(<ut  -  kz  +  Sy  )  ’ 

the  Stokes  vector  components  are  defined  as  follows: 


fv  \ 
*->0 

f  Al + Al  \ 

A, 

Al  ~Al 

— 

x  y 

2  AxAy  cos  8 

UJ 

v2 A.xA.y  sin  8 y 
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where  S=Sy-SK.  It  can  easily  be  checked  that 

Sq  =  >Sj2  +  Si  +  Si ,  (18) 

so  all  the  light  polarization  states  belong  to  the  sphere  surface  of  radius  So  (it  is  called  Poincare 
sphere),  which  is  just  represent  the  intensity  of  the  polarized  light.  The  points  on  the  equator  of 
the  Poincare  sphere  represent  all  the  linear  polarizations,  while  the  other  ones  are  for  elliptical 
polarizations,  which  become  circular  at  the  poles.  According  to  the  classic  definition,  the  circular 
polarization  is  right  handed  if  the  clockwise  time  revolution  of  the  field  vector  takes  place  for  an 
observer  the  photon  arrives  to.  If  this  convention  is  used  then  the  right-handed  circular  state  is 
located  at  the  north  pole  of  the  Poincare  sphere.  For  the  set-up  shown  in  Fig.  5  the  light 


polarization  states  should  belong  to  the  90°  meridian,  and  the  points  where  this  meridian  is 
crossed  with  the  equator  correspond  to  the  linear  polarization  states,  which  are  either  completely 
blocked  or  transmitted  by  the  output  polarizer. 

3.3  Results  of  the  simulations 

3.3.1  Transmission  response  and  dynamics  of  the  polarization  state  in  absence  of 
hydrodynamic  flows 

The  results  of  numerical  simulations  are  shown  in  Fig.  6  and  Fig.  7.  The  Fig.  6  shows  the 
transmission  coefficient  change  under  low  frequency  (t=  0-2  ms)  and  high  frequency  (t=2  -  4 
ms)  driving  with  an  external  voltage  of  50  V  (effective  value)  at  A,=550nm.  The  total  phase 
retardation  change  is  somewhat  less  than  4ti.  The  response  to  the  high  frequency  pulse  is  slower 
compared  to  the  low  frequency  one  for  two  reasons.  Firstly,  this  is  due  to  a  lower  value  of  the 
high  frequency  anisotropy  (sa  =  -4)  compared  to  the  low  frequency  value  (sa  =  +5).  Secondly,  the 
director  field  distribution  at  the  moment  the  high  frequency  pulse  applied  is  very  close  to  the 
homeotropic  state  and,  thus,  very  low  torque  appears  at  the  beginning  of  high  frequency  pulse,  so 
significant  delay  of  500  microseconds  is  well  visible  during  interval  from  t=2  to  t=2.5  ms. 


Fig.  6.  Dynamics  of  optical  transmission  response  in  case  of  dual-frequency  driving,  when 
the  hydrodynamic  coupling  is  neglected 

It  is  very  important  that  the  director  reorientation  is  stable.  The  molecular  reorientation 
under  the  electric  field  appears  just  in  the  xz  plane  that  results  in  the  light  polarization  states 


represented  by  points  on  the  Poincare  sphere  shown  in  Fig.  7.  The  points  representing  the  light 
polarization  state  belong  to  the  90°  meridian  (Si=0).  The  whole  27i-tum  over  the  meridian 
corresponds  to  2n  optical  phase  retardation.  In  our  particular  case  we  have  almost  two  complete 
turns.  The  reorientation  remains  stable  even  if  we  introduce  a  “technological”  azimuthal  offset  in 
easy  axes  orientations  of  about  few  degrees  at  the  two  LC  layer  boundaries. 


S3 


Fig.  7.  Light  polarization  states  during  high  frequency  driving  (t=  2  -  4  ms  in  Fig.  6),  when 
the  hydrodynamic  coupling  is  neglected 


3.3.2  Back  flow  effect  to  the  transmission  response  and  dynamics  of  the  polarization 
state 

The  hydrodynamic  coupling  dramatically  influences  the  stability  of  the  phase  retardation 
switching  under  dual-frequency  driving.  This  is  demonstrated  by  simulation  results  shown  in 
Fig.  8.  The  influence  of  the  hydrodynamics  is  especially  visible  during  the  high  frequency 
driving.  In  this  case  the  transmission  response  does  not  show  the  characteristic  oscillations  of  the 
transmitted  light  intensity  similar  to  ones  shown  in  Fig.  6,  when  no  hydrodynamic  flow  appears 
in  the  system.  The  light  polarization  states  before  the  output  polarizer  also  show  unusual 
behavior,  Fig.  9.  Now  the  points  on  the  Poincare  sphere  do  not  belong  just  to  the  90°  meridian. 


single  package  dual  frequency  driving;  (ULF=UHF=50V); 
(a2,  a4,  a5)=(-0.5,  0.25,  0.25)  Pa  s;  AsLF=+5;  A8HF=-4; 


Fig.  8.  Dynamics  of  optical  transmission  response  in  case  of  dual-frequency  driving,  when 
the  hydrodynamic  coupling  is  taken  into  account. 


Fig.  9.  Light  polarization  states  during  high  frequency  driving  (t=  2  -  4  ms  in  figure  8),  when 
the  hydrodynamic  coupling  is  taken  into  account. 


Only  at  start  of  the  high  frequency  pulse  the  polarization  state  points  revolve  around  Si  axis 
by  movement  along  the  90°  meridian  from  point  at  (Si,  S2,  Sj)=(0,  1,  0)  to  the  “north”  pole  (Si, 
S2,  S:,)=(0,  0,  1)  and  slightly  more.  Further  polarization  states  become  characterizing  by  non  zero 
value  of  the  Si  parameter,  and  they  significantly  deviate  from  the  90°  meredian.  Also  no 
complete  turn  of  the  polarization  states  on  the  Poincare  sphere  exists,  and  thus,  even  271- 
modulation  of  the  phase  retardation  is  not  realized. 

The  results  in  Fig.  8  and  Fig.  9  are  found  for  012/04,5=2.  We  also  found  that  at  higher  ratios 
0(2/014,5  the  numerical  solution  becomes  extremely  unstable,  while  decreasing  this  ratio  results  in 
stabilizing  the  optical  transmission  response  and  corresponding  revolution  of  the  polarization 
states,  Fig.  10  and  Fig.  1 1.  As  we  see  at  0(2/04, 5=5/3  optical  retardation  modulation  close  to  2n  is 
already  realized.  Moreover,  the  decreasing  high  frequency  voltage  (30V)  allows  almost  3  n 
modulation.  All  this  is  well  visible  on  the  Poincare  sphere  in  Fig.  11,  where  one  can  see  the 
revolution  of  the  polarization  states  around  the  Si  axis  along  the  curves  very  close  to  the  90° 
meridian. 


single  package  dual  frequency  driving;  (ULF=UHF=50V); 
(a2,  a4,  a5)=(-0.5,  0.3,  0.3)  Pa  S  ;  AsLF=+5;  AsHF=-4; 


t,  ms 


Fig.  10.  Dynamics  of  optical  transmission  response  in  case  of  dual-frequency  driving.  The 
hydrodynamic  coupling  is  taken  into  account.  On  the  interval  of  high  frequency 
driving  the  response  is  shown  for  two  values  of  the  effective  voltage. 
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Fig.  11.  Light  polarization  states  during  high  frequency  driving  t=  2  -  4  ms  (see  figure  10)  at 
UHf=30V. 


UHF-driving;  (ULF=UHF=50V); 

(a2,  a4,  a5)=(-0.5,  0.25,  0.25)  Pa  S 


Fig.  12.  Hydrodynamic  flow  velocities  (the  x-  and  y-  components)  during  high  frequency 
driving  t=  2  -  4  ms  (see  figure  8)  at  Uhf=50V. 


The  evolution  of  the  back-flow  velocity  shown  in  Fig.  12  allows  complete  understanding  the 
restrictions  for  fast  phase  modulation  coming  from  hydrodynamics.  Indeed,  in  presence  of 
thermal  fluctuations  or  small  azimuthal  off-set  (it  can  be  much  less  than  one  degree)  of  the  easy 
axes  orientations  at  the  two  boundaries  the  flow  is  not  more  characterized  by  single  Vx 
component  of  the  velocity.  The  Vy  component  characterized  by  very  high  values  (~5  mm/s  ) 
appears.  This  flow  in  the  y-  direction  is  the  main  reason  for  disturbing  the  director  reorientation 
making  it  out  of  the  xz-plane  and  providing  the  phase  retardation  behavior  discussed  above. 

It  is  very  interesting  that  there  is  a  small  thickness  range  close  to  each  of  the  layer  boundaries 
where  the  y-component  of  the  flow  velocity  remains  sufficiently  low.  As  one  can  see  from  Fig. 
12  this  thickness  range  is  of  about  20%  of  the  total  LC  layer  thickness.  Thus  for  the  fast  and 
stable  retardation  switching  namely  this  thickness  range  should  be  used.  It  allows  us  to  formulate 
a  simple  practical  rule  for  an  estimation  of  the  optimal  LC  layer  thickness  for  the  fast  271- 
modulation  as: 


d  = 


A 


(19A.17) 


0.2An 

Thus,  for  example,  at  A»=0.25  the  optimal  cell  gap  for  T=633nm  should  be  of  about  12.5  pm. 


4.  Conclusions 

Using  numerical  simulations  we  have  established  the  main  mechanism  providing  restrictions 
for  fast  retardation  switching  that  is  due  to  by  the  hydrodynamic  destabilization  of  the  plane 
director  field  distribution.  In  particular,  using  Poincare  sphere  presentation  we  have  studied  the 
dynamics  of  the  light  polarization  states  in  absence  and  presence  of  the  hydrodynamic  coupling. 
It  is  also  shown  that  the  hydrodynamic  destabilization  can  be  suppressed  by  increasing  the  ratio 
of  Leslie  coefficients'  1 012/04,5 1.  Also  we  have  found  that  close  to  the  LC  layer  boundaries  there  is 
space  region  of  20%  of  the  LC  layer  thickness  d,  where  the  hydrodynamic  destabilization  is 
significantly  lower  compared  to  the  rest  region  in  bulk  of  the  LC  cell  gap.  Thus  namely  this 
effective  thickness  should  be  used  for  the  fast  modulation  of  the  optical  retardation. 
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Figure  captions 


Fig.  1.  Electro-optical  response  measured  for  the  0-aligned  AP-99-1  layers  in  RP45  and  RPo 
geometries  at  Ulf  =  50V  and  T  =  22  C):  a)  d  =  7.6  pm,  curve  1-4  are  for  Uhf  =22,  30,  40  and  0 
V,  respectively;  b)  d  =10  pm,  curve  1-3  are  for  Uhf  =  30,  38  and  51V,  respectively.  The  lower 
trace  represents  the  applied  voltage  waveform. 

Fig.  2.  Thickness  dependence  of  high  frequency  (35  kHz)  voltage  Uhf  at  which  a  dynamic  phase 
retardation  change  of  n  (3)  and  2k  (2)  is  still  achieved  (100%  light  modulation  efficiency  is 
observed).  The  curve  1  shows  dependence  of  the  threshold  voltage  of  dynamic  phase  modulation 
breakdown.  The  data  are  for  the  0-aligned  AP-99-1  layers  at  T  =  22  C. 

Fig.  3.  A  sketch  of  the  double-layer  liquid  crystal  cell  and  its  connection  to  the  net  source.  The 
arrows  show  the  directions  of  rubbing. 

Fig.  4.  Electro-optical  response  of  the  0-aligned  AP-99-1  layers  in  a  double  layer  (6.8  and  7.3 

O 

pm)  cell  measured  in  RP45  and  RPo  geometries  at  Ulf  =  50V  (T  =  22  C):  curve  1-4  are  for  Uhf 
=0,  22,  28  and  32  V,  respectively. 

Fig.  5.  Scheme  of  simulated  electro-optical  set-up. 

Fig.  6.  Dynamics  of  optical  transmission  response  in  case  of  dual-frequency  driving,  when  the 
hydrodynamic  coupling  is  neglected 

Fig.  7.  Light  polarization  states  during  high  frequency  driving  (t=  2  -  4  ms  in  Fig.  6),  when 
hydrodynamic  coupling  is  neglected 

Fig.  8.  Dynamics  of  optical  transmission  response  in  case  of  dual- frequency  driving,  when  the 
hydrodynamic  coupling  is  taken  into  account. 

Fig.  9.  Light  polarization  states  during  high  frequency  driving  (t=  2  -  4  ms  in  figure  8),  when  the 
hydrodynamic  coupling  is  taken  into  account. 


Fig.  10.  Dynamics  of  optical  transmission  response  in  case  of  dual- frequency  driving.  The 
hydrodynamic  coupling  is  taken  into  account.  On  the  interval  of  high  frequency  driving  the 
response  is  shown  for  two  values  of  the  effective  voltage. 


Fig.  11.  Light  polarization  states  during  high  frequency  driving  t=  2  -  4  ms  (see  figure  10)  at 
Uhf=30V. 


Fig.  12.  Hydrodynamic  flow  velocities  (the  x-  andy-  components)  during  high  frequency  driving 
t=  2  -  4  ms  (see  figure  8)  at  Uhf=50V. 


